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o 

Abstract. We give a direct and simple proof of Touchard's continued fraction, pro- 
vide an extension of it, and transform it into similar expansions related to Motzkin 
f T | , and Schroder numbers. Another proof is then given that uses only induction. We use 

■ this machinery on two examples that appear in recent papers of Josuat- Verges; with 

, an additional parameter, these two can be treated simultaneously. 

o 

U; 1. Introduction 

Touchard [8] studied a certain function F that he developed as a continued fraction 
(equation (20), loc. cit.): 



(N 
O 



l-z(l-q)F(q;z(l-q)) ^ (l-q)z 



, (1 - <? 2 )^ 

00 1 1 

^ ■ He (equation (32), loc. cit.) also found the identity 



1 



v ' k>0 



l- Z (l-q)F(q;z(l-q)) l-zC( 
where C (z) is the generating function of the Catalan numbers: 



£ : c{z) = — 2z — 

n>0 



>^ 1 (2n\ 
^ n + 1 1 n J 



Thanks to the substitution z = nj^p > which implies C(z) = 1+v, Touchard's continued 
fraction takes the following form (a "T-fraction"): 

-=J2(-l) k q^v k . (1.1) 



(1 — q)v 

1 + v — 

(l-g> 
1 + v 

Touchard's combinatorial approach (counting chord diagrams) is ingenious but at the 
same time quite involved. The aim of the present note is to provide a simple and direct 
proof of i jTTjl . 
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This project was started in June 2009 and was planned to be a cooperation between 
the present author and Philippe Flajolet. Unfortunately, this never materialised, and 
it was continued in 2011 with more continued fractions, taken from papers by Josuat- 
Verges [Hill]. 

There is apparently a rich combinatorial world behind these continued fractions. 
However, here, the emphasis is to leave the combinatorics completely out and do 
"purely manipulative" proofs. 

Perhaps the simplest proof is by induction, if one "knows" already the result for the 
continued fraction with the first i lines deleted. This point of view is elaborated on 
two examples. 

We did not go for an exhaustive list of continued fractions of the Touchard type, 
in order to keep this paper short and crisp. However, we are confident that the ap- 
proaches) taken here would also work in other similar instances. 

This paper is not overlap-free from other papers, in particular [31 0], but we are 
confident that it contains enough original material. 

2. The proof 
We introduce a parameter t and consider 

*w- 1 



(l-qt)v 1 + v - (1 - qt)vF(qt) ' 

(1 - qH)v 
1 + v 

Setting F(t) = A(t)/B(t), this leads to A(t) = B(qt) and 

B{t) = (1 + v)B(qt) - (1 - qt)vB(qH). 

We make the ansatz 

B(t) = ^ 

n>0 

and get upon comparing coefficients 

/i i \ n 2n , 2n— 1 V Q a n-l 

a n = (1 + v)q a n - vq a n + vq a n _i = T — -r. 

(1 — q n )(l — vq n ) 

This can be iterated, and since ao = 1, 

2 

v n q n 



(q;q)n(vq;q)n 

where we employed the notation (x; q) n = (1 — x)(l — xq) ... (1 — xq n ~ l ). Thus 

^ {q; q)n{vq; q) n 

and Touchard's continued fraction is given by F(l) = B(q)/B(l). 
Our goal will be achieved once we are able to establish the identity 

fk+i\ 



B(q) = £(1)5>1)V 



2 J V K . 



k>0 
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This is the same as 



E 



v q 



E 



n > (q;q)n(v,q) n ^ (?; ?)„(«; g)„ n > 

where we replaced t> by f/g. One of the sums can be evaluated thanks to a formula 
due to Cauchy [H (10.9.2)]: 



E 



o (q;q) n (v;q) n ^ (?5?)n' 



E 



Comparing coefficients of t> ra , we are left to prove that0 



n — 1 
k-1 



E 



Here, we use the notation 



1 k=0 



(?;?)n 



:-i)v 



To establish such an identity is nowadays routine, and the tool for it is the g-version 
of Zeilberger's algorithm [6|. It produces the recursion 

(g n+1 - g)T(n) + {q 2n+3 - q n+2 + q- q n+l + l)T(n + 1) + (g n+2 - l)T(n + 2) = 0, 

for both sides, and together with a few easily checked initial conditions, we are done. 



3. An extension 
Our proof gives more, by considering F(q l ): 

B(q l+1 ) 1 



F(q l 



l + v 



1 — q % )v 



l + v - 



This quotient by itself it not nice, but we have 



B(q 



IT^E-D" 



B(l) 



n>0 



n + i 
i 



( n+1 ) n 



The identity behind this is 



E 



(?)n(*0n (?)n 



Et^EM)"""^ 



n + « 
i 



which, in terms of coefficients, reads as 



n „k 2 +ik 

V- 



n — 1 
fe-1 



E 



i k=0 



(?;?) 



n—k 



k + i 

i 



X J. Cigler, who has seen an early draft of this paper, has shown me a "human" proof. 
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Again, Zeilberger's algorithm proves this readily. Therefore 
B(q i+1 ) I B(q i ) 



F(q l ) 



B(l) I 5(1) 
n + i 



£(-i) B 



n>0 



"7D- 1 )' 



n>0 



n + i — 1 
i - 1 



4. Another proof, starting from the righthand side 



We set 



n>0 



n + i 
i 



qK 2 )v n i> 0, 



so 



S-i := 1 and expand the quantity as a continued fraction of the form 

S-l 



l + v - 



1 + v - 



a 2 v 



It is not hard to show that the sequence 01,02,... is unique, as one can compare 
coefficients in series expansions of both sides, and compute them one by one. 
We claim that Oj = 1 — q % and prove the formula 



1 + v - 



a i+1 v 



Si-l 



1 + v - 



by induction. Since So = 1 — qv + ■ ■ ■ , the value a± is established. And now 

1 Si 



1 , S *+ 1 S- 1 

l + v — a i+1 v s*- 1 

Si 

leads to 

(1 + v)si - a i+1 vs i+1 = Si-i. (4.1) 

(All the steps are reversible, and it might be clearer to start the arguments from this 
recursion.) Comparing coefficients, 

n + i — 1 

i 

n + i — 1 
i-1 



(-1)' 



n + % 
i 



Q 



= (-i) b 



q 



q {2) _(l_g»+l)(_l) 

J 



n + i 
i + 1 



q 



or equivalent ly 



n + i 




n + i — 1 










9 


i 



+ (l-q 



n + i 




n + i — 1 




9 


i-1 



q n , 
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or 



or 



n + i 




n + i — 1 




q n - 




i 


i 


% 



n + i — 1 




n + i — 1 


i 




z - 1 



n + i 


-q l 

i 


n + i — 1 




n + i — 1 


i 
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i - 1 



which is the basic recursion; thus the proof is finished. 
Remark. J. Cigler has kindly pointed out that 



£ 

fc>0 



-i - 1 
k 



1/9 



which are Rogers-Szego polynomial with negative indices. 

Furthermore, for v = —1 (14.11) simplifies, and the recursion can be solved by iteration, 
leading to a product representation. 

5. MOTZKIN NUMBERS 

Cigler [2] recentlj^l found the identity 

I «, . (k+l\ 2k 1 \ / r/\2k+l 



1 - Z 



(l-q)z 2 



k>0 



;i - q 2 ) z 2 



l-z 



where 



M{z) 



l-z- y/l-2z- 3; 
2? 



is the generating function of the Motzkin numbers. 

Using the substitution z = 1+ " +v 2 , we see that, after division of both sides by 
1 +v + v 2 , only even powers of v appear, and with the further substitution v 2 = x, the 
formula transforms into Touchard's continued fraction. 

6. Schroder numbers 
Cigler also found the formula 

1 



1 - z 



(l-q)z 



J>i)V'Vs(*) 



2k+l 



k>0 



1 - Z 



l-g 3 * 



where 



S(z) 



l-z-^l-6z + z 2 
2z 



is the generating function of the Schroder numbers. 



2 In 2009, when a first draft of this paper was sketched. 
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We can also reduce this to an instance of Touchard's continued fraction as follows: 
The substitution y — z/(l — z) 2 transforms the formula into 



1 - 



1 - 



(i - q)y 
(i-q 2 )y 



z )J2(-l) k q^z k S(z) 2k+1 . 



k>0 



But recall that 



(l-q)y l-yC(y) 



(i-q 2 )y 



y vw; fc>0 



so that the claim follows from the equalities 

(1 - z)S(z) ' 



and zS 2 {z) = C{y) - 1. 



1 - yC{y) 

But both are a routine verification that is best done with a computer algebra system. 



Writing 



7. A FORMULA DERIVED BY RlORDAN 
1 



1 - zF(q; z) 
we find, upon comparing coefficients, that 

T n (q)(l - q) n = [z n ] - - 1 Yl i 1 - C ^)' 



k>0 



[v n ](i-v)(i + v) 2n J2 q y* )(-v) 

(k+i\ \ ( In 



q 



k=0 



n — k 



2n 
n — k — 1 



which is a formula derived by Riordan [7]. 

8. Another continued fraction 
We consider now an example from j3l H] : 



1 + v- 



(l-q) 2 v 



1 + v- 



(l-q 2 ) 2 v 



We set 



F(t) 



1 + v- 



[l-qt) 2 v l + v-{l-qt) 2 vF(qt)' 



1 + v- 



(1 - q 2 t) 2 v 
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This continued fractions appears starting from 

1 



Cl-qfzJ 
1 (l~q 2 ) 2 z 



after the usual substitution. Set 

Fit) = W 

{) B(ty 

We get A{t) = B(qt) and 

B(t) = (1 + v)B(qt) - (1 - qt) 2 vB{q 2 t). 

Reading off coefficients of t n leads to a recursion of second order, which is not nice. 
But if we define 

«*) - . 

(*; ?)oo 

then 

(1 - t)/3(t) = (1 + - (1 - gt)t;/3(g 2 t), 

and 

Pn ~ Pn-1 = (1 + " + g 2 "" Vn-1, 

or 

1 + vq 2n ^ (-vq; q 2 ) n 

Pn — Pn-r 



(1 - q n )(l -vq n ) (g;g) n (wg;g) ? 
since /3 = 1. Hence 

(g; q) n {vq; q) n 

and 



(l-t)jS(t) 

We need F(l). The interpretation of (1 — t)(3(t) IS cLS cL limit. We use the Heine 
transform, as was done in [3j 0] as well: 



^ (q;q)n(vq;q)n 



and therefore 



iv/ug; g) (iy^g; g) n (*; g) 

pit) = -. r r - 51 r > - r— (— lWVq) n 

{vq; q)oo{t; g)oo ^ (g; q) n {iy/vqt; q) n 

and 

(1 - t)p(t) = ( i v / ^^ q)ooi.-i\f v li g)oo (iy^g; g)n(t; g)n . ,„ 
(«g; g)oo(^g; g)oo f^{q;q) n {i^vqt;q) n 



Now we can plug in £ = 1: 



g)oo(g;g) 
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Furthermore 



or 



Hence 



(w, 2)00(2; g)oo ^ (q;q)n(i^vgq;q)n 
F(l) = V — (-iJvqT 

n,m>0 

= i m+n (^/vq) m+n q nm (-l) n 

n,m>0 

= £ (-i)"( ug )V (2JV ~ n) (-i) n - 



0<ra<2iV 



So 



[^]F(l) = £ (^AT+n^nW 
0<n<27V 

_ ^ ^_ 1 ^n^Af+2(n+7V)Af-(r i +Af) 2 
-AT<n<Af 



-Af<n<Af 

To find the coefficients (in the variable z) of the original continued fraction, we compute 
[^](1 + v)F(l) = [v N ]{l + v) 2N (l - v)F(l) 

N 

= + v) 2N [v k ](l - v)F(l) 

k=0 

N 



k=0 L v 7 v 7 J -k<j<k 



This translation from the coefficients in the variable v to z works always in exactly the 
same way, so we will not state the z-formula for the further examples. 

As a bonus, we can also consider 

1 



(1 - q^v 

1 + v 



1 + v 

which can be computed via 

B(q i+1 ) B(q t+1 ) I B{q i ) 



B(qi) ~ 5(1) / 5(1)' 
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Now 
B(l) 



{-vq;q 2 )oo 
(vq;q)oc 



= i (iy^g; q)n(g l ; q)n | 

{\^/vq\ q)i ^ (q; g)„(i v /ugg i ; g)„ 



-iv^) r 



n>0 



(g; g)n(i-v/^9; 9) 



n+i 



1 »,^n v 1 



(9; ^ (iv 7 ^^; s) 



(9; ?) 



i-l 



i + m — 1 

m 



(9; 9) 



i-l 



E 

n,m>0 

E<>" E 

iV>0 0<n<2iV 



(g^ 1 ;g),_ 1 (i v ^O m (-i v ^) n 



i + 2N — n — 1 
i - 1 



— r 



(9; q)i- 



v N g N(N+l) 



E 



iV>0 



N>0 



-N<n<N 



-N<n<N 

'i + N-n-1 
i-l 



'i + N-n-1 
i - 1 



(q n+1 ;q)i-iq N+2nN - n2 (-l) N+n 

(? B+JV+1 ;?)i-i9- na (-i) T 

g" n2 (-l) n 



i + iV + n- 1 
% - 1 



and the continued fraction has been evaluated as a quotient of two series of this type. 

9. Independent proof by induction 

As before we can turn this extra information into an extremely elementary proof by 
induction: We set 



N N(N+1) 



> r "(I ^2 

N>0 ~N<n<N 
SO 



~i + N -n 




r i + N + n 


i 


q 


i 



q- n (-1)", 



s_i := 1, and expand the quantity as a continued fraction of the form 

s-i 



l + v - 



CL\V 



l + v - 



a 2 v 



We claim that at = (1 — q l ) and prove the formula 

1 Si 



1 + v - 



a i+1 v 



Si-l 



l + v - 



a i+2 v 
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by induction. Since l/s Q = 1 — q(q — 2)v + • • ■ , the value a± is established. And now 

1 Si 



1 i Si + 1 S- i 

1 + v — a i+ iv 



leads to 
Set 



Si,N ■— Q / J 



-N<n<N 



i ~ a>i+lVSi- 


hi 


= Si_i. 


~i + N — ri 




"z + iV + n 


i 


i 


i 



(9.1) 



<T" (-1)' 



Then we must show the above recursion, which, in terms of the coefficients, reads as 

Si,N + Si,7V-l — (1 — q %+ ) 2 Si+l,JV-l — Si-i s N = 0. 

Manuel Kauers (Rise, Linz) kindly provided an automatic proof for this, using the 
package HolonomicFu notions by Christoph Koutschan [5]. 

Remark. J. Cigler has pointed out that (19. ip becomes nicer when written with nega- 
tive indices: Set Si := s_j_i, then 



S i = (l + v)S i - 1 -(l-q 1 - i ) 2 vS, 



i-2i 



and 



S = 


1, 




% 




i 


_N-n 


9 





iV=0 -N<n<N 

For g = 1, this is just the binomial theorem. 

10. Another continued fraction of the Schroder type 
In this section we investigate the following continued fraction: (cf. [3l (18)]) 

1 



1 + v- 



(l-g)(l-g> 



1 + v - 



(1 - q 2 )(l - q 3 )v 



Set 



F(t) 



1 + v 



(1 - qt)(l - qH)v 1 + v - (1 - qt)(l - q 2 t)vF(qt) ' 



(1 - q 2 t){l - qH)v 



1 + v 



With the usual F(t) = A(t)/B(t), we find A(t) = B(qt) and 

B(t) = (1 + v)B(qt) - (1 - qt)(l - q 2 t)vB(q 2 t). 
We set p(t) = fl(t)/(gi; ?)«,: 

(1 - gt)/3(t) = (1 + - (1 - qt)v(3(q 2 t). 

Hence 

/3 n - qPn-1 = (1 + %"/?n " V"/?n + V^Pn-l, 
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or 



q(l+vq 2n ~ 2 ) q n (-v;q% 

-Pn-l — 



(l-<r)(l-v<r) yn (q;q)n{vq;q) n 



So 



Hence 



= (-iy^; g)oo(iVt;gt; g)^ (iy^g; g) n {qt; g) n , . 
(vq;q)oc{qt;q)oc ^{q;q)n(SVvqt;q) n 



m _ (-«;?% 



(w, q)oo(q; q)oo l + iv^ 



and 



o/ x = (-iy^;g)oo(ivW;g)oo (V^g;g)n(g 2 ;g)n , . 
(vq;q)oc(q 2 ;q)oo ^ (q;q)n(iVvq 2 ;q)n 

(-v,q 2 U v (l-g w+1 ) 



(vg; g)oo(g; g)co ^ (l - iv^g n+1 



Therefore 

5(g) (g 2 ;g)oo/?(g) 



F(l) 



l-q Z-.(l-ufcq n+l Y V ^ 



f+n v Nq(n+l)(2N-n) 



5(1) (g;g)oo/3(l) 

(i-q 

(1 - V^g n+1 ) 
= E (! - ^ +1 )(iv^g n+1 ) m (-iv^) n 

^ n,m>0 
^ Af>0 0<n<2AT 

+ 1 ^ (1 _ g n+l^_ 1 ^+n+l u iV+l 5 (n+l)(2JV+l- 

^ ^ JV>0 0<n<2iV+l 

= 1 ^ v N q NiN+1) ^ i 1 ~ q n+N+1 )(-l) n q~ n{n+1) 

^ 1 N>0 —N<n<N 

-^Lj^v 2 Yl (i-? n+JV )(-i)v- 2 . 



n>0 



iV>0 -N<n<N 

One could rearrange that, but it would not get any better. 
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11. Generalisation of the last two continued fractions 
Here we consider 



(l-g)(l-g> 

^ (l-g 2 )(l-g d +> 
1 + 1> 

with d a nonnegative integer. So for d — 1 and d = 2 we get the continued fractions 
from the previous sections. As usual 

m = 77 FS^S ' 

1 + v 



(1 - qH)(l - q d +H)v 
1 + v 



so 

F(t ) = I = 

W 1 + (1 - qt)(l - q d t)vF(qt) B{q) 

with 

B(t) = (1 + v)B(qt) - (1 - gi)(l - q d t)vB{q 2 t). 
Set p(t)=B(t)/(q d -H; q)^, then 

(1 - g^)/3(t) = (1 + W )/3(gt) - (1 - qt)v(3(q 2 t). 

Comparing coefficients, 

Pn ~ q^Pn-i = (1 + - vq 2n $ n + V^X-i, 



or 



g^l + vg") g( d - 1 ) n (-ug 2 - d ; g 2 ) r 

Pn-l — 



So 



(1 -g n )(l -vq n ) (q;q)n(vq;q) n 

(?;?)n(«?;?)n (9 j 



n>0 



N; g)oo(g d - 1 *; g)oo ^ (2; V^g d/2 *; g)n g 

Let us exclude the case d = 1, so that we don't have to worry about taking a limit. 
5(1) = {-vq 2 - d ; q 2 )oo v- (g d " 1 ;g)n ( _. i- d/2) n 

^ = (-V- d ;g 2 )oo \- (iy^g d/2 ;g)n(g d ;g)n , ■ /r. n i-d/2 )n 

PW (vq; g)oo(g rf ; g)oo(iv / ^g 1 - <i/2 ; g)<*-i ^ (g; g)„(ivW +c!/2 ; g)n + i { v q ) 

(-vq 2 - d ;q 2 ) 00 



(w, g)oo(g d ; g)oo(ivW~ d/2 ; g)d-i 

(1 - VUg d / 2 )(g d ; g) 



V 1 -iy/uq- )yq ,_qjn r~ l-d/2 )n 

(g; g) n (l - iv^g 1+d / 2+n )(l - i^q d l 2+n ) 1 V q ' ' 



X 

n>0 
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So 



F(l) = ^ 
U (1 — q d ~ l )(3{l) 



V (g d ;g)„(l -iy^g d/2 ) , /T, n 1 - d / 2 ) n 

^(g;g) n (l-i v ^g 1 + d / 2 +")(l-i>^ /2+n ) 

(_s /r l-d/2\n 

, y (q d ;qUl-q n+1 ) . r-i-d^n 



n>0 



£ ?aSf^' i ^ /2+ "» m (- i ^ 1 " J/2 »" 

n>l,m>0 ^ v U\wjn i 

+ E (9 In )n( \7 (ivW +d/2+ "H-ivW- d/2 y 



{q d ]q)n{ l) N+n ^ n (2N-n)+Nd~nd 

(1 - q)(q;q) n -i 



JV>0 l<ri<2Af 

, JV V"^ (g^; g)n(l - g n+1 )(-l) jV+w (i +n )(2N-n)+Nd-nd 



JV>0 0<n<2JV ^ 9){Q',Q)n 
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